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Abstract 

The time evolution of 'fiducial' radio sources derived from fitting 
the dynamical model of Kaiser et al. (1997) is compared with the 
observational data for the 'clan' sources found in the sample of gi- 
ant and normal-size FRII-type sources published Paper I (Machalski 
et al. 2004). Each 'clan' comprises 3, 4 or 5 sample sources having 
similar values of the two basic physical parameters: the jet power Qq 
and central density of the galaxy nucleus po (determined in Paper I) 
but different ages, radio luminosities and axial ratios. These sources 
are considered as the 'same' source observed at different epochs of 
its lifetime and used to fit the evolutionary luminosity-size (P — D) 
and energy density-total energy (u c — E tot ) tracks derived from the 
model for a 'fiducial' source with Qq and po equal to the means of rel- 
evant values obtained for the 'clan' members, as well as to constrain 
the evolutionary model of the source dynamics used. In the result we 
find that (i) The best fit is achieved when the Kaiser et al.'s model 
is modified by allowing an evolution of the sources' cocoon axial ratio 
with time as suggested by Blundell et al. (1999). (ii) A slow acceler- 
ation of the average expansion speed of the cocoon along the jet axis 
is suggested by the 'clan' sources. We argue that this acceleration, 
although minor, may be real and some supporting arguments come 
from the well known hydro dynamical considerations. 

Key words: galaxies: active - galaxies: evolution - galaxies: 
kinematics and dynamics 

1 Introduction 

There are several analytical models of the dynamical evolution of extended 
powerful radio sources whose prototype is the radio galaxy CygA (3C405). 
All of them are based on the hydrodynamical expansion of their lobes (or 
cocoon) caused by the interaction of light supersonic jets with the ambient 
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material swept-up within the bow-shock arising behind the jet head, and 
satisfactory describe the sources' dynamics. However, only two of them, the 
most sophisticated ones published by Kaiser, Dennett-Thorpe and Alexander 
(1997) [hereafter KDA] and Blunder!, Rawlings and Willott (1999) [hereafter 
BRW], deal with the sources' energetics and the radio luminosity evolution 
with time. These two models differ in one assumption: the KDA model sus- 
tain the Falle's (1991) requirement of a homologous (self-similar) bow-shock 
expansion, while in the BRW one the jet-head pressure does not scale with 
time in the same way as cocoon pressure, hence it does not expand self- 
similarly. As the result, the predicted decrease of the radio luminosity of 
matured sources is faster in the BRW model than in the KDA one. Besides, 
the self-similar bow-shock expansion would be possible only if the internal 
source (cocoon) pressure was always much higher than the external envi- 
ronment pressure, what has been questioned by X-ray observations of the 
medium surrounding selected FRII-type galaxies and quasars (e.g. Hardcas- 
tle and Worrall 2000). Therefore giant radio sources with their largest linear 
sizes and low luminosity are excellent objects for an observational constraint 
of the both models predictions. 

In Machalski et al. (2004) [hereafter referred to as Paper I] we presented 
the sample of 18 giant lobe-dominated radio sources (with projected linear 
size D > 1 Mpc if H =50 kms _1 Mpc _1 and f2 M =l) and 54 normal-size 
sources selected to fulfill the following criteria: (i) the sources have the FRII- 
type (Fanaroff and Riley 1974) morphology, (ii) the existing radio maps allow 
a suitable determination of their lateral size transversal to the source's axis 
which is necessary to specify its (so called) 'axial ratio' and volume, and (iii) 
their spectral age or the expansion speed, reliably determined with the same 
model of the energy losses, are available from the literature. 

In Paper I, applying the KDA model chosen for its simplicity in com- 
parison to the BRW, for each member of the sample we have determined 
their basic physical parameters, i.e. the jet power Q , the central density of 
the galaxy nucleus po, the energy density and pressure in the lobes/cocoon 
(u c and p c ), and the total energy of the source E tot . All these parameters 
are derived from fitting the model free parameters to the source redshift, 
monochromatic radio luminosity, projected size, and axial ratio, i.e. its ob- 
served parameters. Next, these model and observed parameters have been 
used to search for conditions or circumstances under which extragalactic ra- 
dio sources can reach their extremal linear extent even of a few Mpc, to check 
the values of the internal pressure in their lobes and to compare them to the 
suggested external pressure in the sources' environment, etc. 

As a result, we argued that (i) the giant-size sources do not form a unique 
class of extragalactic objects; they are old sources with the jet power Q > 
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3 • 10 37 W evolved in a relatively low-density environment with p < 10 
kgm -3 , (ii) an apparent increase of the lowest internal pressure value with 
redshift, found in the largest sources, is obscured by the intrinsic dependence 
of their size on age and the age on redshift, (iii) the ratio of the jet-head 
pressure expanding the source cocoon in the jet direction and the cocoon 
lateral pressure (ph/p c ) governing its axial ratio is dependent on the source's 
age (as well as on the jet power) which violates the KDA model assumption 
of self-similar expansion of the cocoon. We showed that a departure from 
self-similarity for large and old sources is justified by the intrinsic correlation 
between their observed axial ratios and ages. This forms an independent 
observational support for the competitive BRW model. 

Moreover, we have realized that for a number of the sample sources the 
derived values of their physical parameters Qo and po were very close, while 
their ages, luminosities, and linear sizes were significantly different. This 
immediately raises a question whether a proximity of the above two funda- 
mental model parameters appears by chance or it has a real physical mean- 
ing. If the latter was true, observed luminosity and size of these sources 
should follow more or less the model time evolution predictions. To check 
this we assume that they may be considered as 'the same' source observed 
at different epochs of its lifetime. Hereafter we call these sources as a 'clan' 
although the members of each 'clan' have rather different redshifts. The 
above may arise a doubt about the density p being constant in the KDA 
model during the source's lifetime at a given redshift. However, there are 
no evidences for an evolution of galaxies' density and density profile during 
the life of a radio source , thus the neutral assumption that such parameters 
are independent of redshift is commonly accepted (cf. Gopal-Krishna and 
Wiita 2001). Nevertheless, the galaxy halo can presumably be scaled with 
the mean matter-density evolution of the Universe proportional to (1 + z) 3 . 
To check how this effect can influence the source's parameters derived from 
the KDA model, we have selected two giants with the largest and smallest 
redshift in the sample, i.e. MRC0437-244 and PKS0319-454 with redshifts 
of 0.84 and 0.0633, respectively (cf. Table 1 of Paper I). Using the cosmolog- 
ical constants adopted in this paper, both sources have been shifted in the 
look-back time by their given ages (19 and 180 Myr, respectively; Table 2 
of Paper I). Then, we have computed the redshift at the birth-time of these 
giants, and recalculate their model parameters for such higher redshifts. The 
resultant parameters differ only very little from the original ones (by less 
than 1%) justifying the assumption that the density profile of a galaxy does 
not evolve during rather short lifetime of radio sources in comparison to a 
time-scale of the matter large structures' evolution. 

Hereafter we show that these 'clans' appear crucial for the observational 
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constraint of the KDA model predictions of the source's time evolution pre- 
sented in this paper. In Section 2 the members of the three 'clans' are spec- 
ified, and the evolutionary tracks of their model representation on the radio 
luminosity-linear size (P — D) and energy density-total energy (u c — E tot ) 
planes are calculated and compared the observational data. The results ob- 
tained, especially the necessity of a departure from self-similarity of the lobe 
(cocoon) expansion in time, are discussed in Section 3. The conclusions are 
given in Section 4. 

2 Evolutionary Tracks of Sources 

All dynamical models enable one to calculate evolutionary tracks of time- 
dependent secondary parameters of radio sources if their primary (time- 
independent) parameters are given. In the papers of KDA and BRW the 
tracks of radio luminosity P versus linear size D were derived for imaginary 
sources with assumed values of Q , p , a , /3, and z. 

In our approach we are able to calculate such evolutionary tracks for 
actual sources. The sample of radio sources used in this analysis (cf. Paper 
I) consists of the four sets of sources: giants (with D >1 Mpc), normal-size 
high-redshift and low-redshift (with z >0.5 and z <0.5, respectively), and 
low-luminosity (with P 14 < 10 244 WHz _1 sr _1 ). In Section 1 the 'clans' of 
sources have been introduced, i.e. the sample sources with very close values 
of the two fundamental parameters in the KDA model, the jet power Q 
and central core density po, and evidently different ages, luminosities, and 
axial ratios. Since this dynamical model assumes constant jet power during 
a source lifetime, and the nucleus density p is a priori constant, members 
of such a clan can be considered as 'the same' source observed at a number 
of different epochs throughout its life. The observed parameters of these 
members can verify predictions of the model. 

2.1 Clans 

We found 6 clans consisting of 3, 4 or 5 sample sources fulfilling the adopted 
selection criterion, namely that the fitted values of Qo and po do not dif- 
fer from their mean value by more than 30 %, i.e. the standard deviation 
from the mean is less than 0.114 in the logarithmic scale. Also the redshift 
of members should be comparable, however in our limited sample we have 
accepted redshift ratios up to about 3. Leaving out three-member clans, the 
remaining three of six clans are marked in Fig. la and b, and their members 
are listed in Table 1. 
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Figure 1: a) Plots of the jet power Q against central density of the core 
p - The giants are indicated by crosses, high-redshift sources - diamonds, 
low-redshift sources - open circles, and low-luminosity sources - stars. The 
'clans' of a few sources with comparable values of Qq and po are marked by 
large circles with the radius being equal to the standard deviation of logQo 
and logp values for given 'clan' members; b) The same 'clans' are localized 
on the u c - Qo plane by time axes. 
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Table 1: Members of the clans 



Clanl 



Clan3 Clan4 



(lg<2o[W]> 
(lgPo[kg/m 3 ]) 



39.26±0.08 

-23.24±0.04 

0.72±0.26 



38.84±0.03 38.21±0.10 
-23.35±0.05 -23.22±0.07 



0.43±0.15 0.20±0.09 
3C166 3C319 



clan members 



3C263.1 

3C289 

3C411 

3C272 

3C274.1 



3C334 1025-229 
1012+488 0136 +396 
0821+695 3C326 



The mean values of logQo, logpo, and redshift in each clan are given with 
the standard deviations of their individual values from the mean. The mem- 
bers are ordered according to their increasing age. In all three clans that 
increase of age is accompanied with an increase of the size and total energy 
of the members, altogether with a decrease of their luminosity and energy 
density. This qualitative observational behaviour, concordant with the model 
expectations, is quantitatively presented in Section 2.2. 



In order to check whether observed parameters of the members of a given clan 
are consistent with the time evolution derived from the KDA model, a fiducial 
source has been created for each clan. This fiducial source has Qo and po 
equal to the mean value of these parameters in a given clan and the evolving 



axial ratio of the cocoon AR(t, Q ) oc £0.23±o.o3g0.i2±o.o2 ^ ^ in Paper ^ 



For each of the clans, a time evolution of its size D(t), luminosity P(t), energy 
density u c (t), and total energy E tot (t) have been calculated. The predicted 
D(t) tracks for the three clans are shown in Fig. 2a with the dotted lines. 
The solid lines connect members of each clan drawn with different symbols. 
The data points are shown with error bars (the errors in size are usually 
smaller than the vertical size of symbols). The dashed line indicates D(t) for 
the constant advance speed of 0.1c. The higher slopes of the modelled tracks 
with respect to this constant speed seem to indicate a slowly accelerating 
expansion of the cocoon in the clans investigated, however the effect is of low 
significance because of uncertainties of the age. More extensive discussion of 
the above effect is given in Section 2.3. 

The time evolution of the total energy predicted with the KDA model for 
the three clans is shown in Figure 2b with the dotted lines. The short solid 



2.2 Time Evolution of the Clans 
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Figure 2: a) Time evolution of the size of the fiducial source for each of 
the three observed clans (dotted lines). The dashed straight line indicates a 
constant advance speed of 0.1c; b) Evolution of the total energy of a fiducial 
source for each of the clans. The dotted lines show the model prediction, the 
shorter solid lines - the best power-law fit to the data points. 
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lines indicate the best power-law fit to the data points. A correspondence 
between the model E tot (t) and the data points is satisfactory. 

2.3 Tracks logP-log£> 

The evolutionary tracks of the fiducial source for the three clans are shown 
in Figure 3. The members of separate clans are marked by different symbols 
as in Figure 2a and 2b. The actual age of each source is indicated by the 
number in a vicinity of each symbol. The dashed curves show the tracks 
calculated from the original KDA model, i.e. with a constant AR taken as 
the mean of axial ratios in a given clan. It is clearly seen that these fits are 
unsatisfactory. Therefore, we make the second run of the fits with the KDA 
model but modified by substitution of Rt in the formula of Kaiser (2000) 

V hc = (2.14 - 0.52^)^ 04 -°- 25/3 , 

used in Paper I, by our empirical formula for the age-dependent -Rt(^) = 
AR/2 oc t - 24 (cf. Paper I). Hence 

V hc = 0.0056 1 04 , 

where t is expressed in [yr] . This is consistent with the formalism used in the 
BRW model in which the injection spectral index is governed by the breaks 
in the energy distribution of the particles injected into the cocoon from the 
hotspot instead of the constant index of p = 2.14 in the KDA model. The 
modified tracks, marked with the solid curves, show that the evolving AR 
much better fits the observed changes of P and D. The markers of the same 
age on these tracks are connected with a dotted line. Further discussion of 
the best fit tracks is given in Section 3.1. 

2.4 Tracks log(w cq )-log(£ tot ) 

The model also allows one to predict the evolution of a source on the energy 
density-total energy plane. These tracks are shown in Figure 4 with dashed 
curves calculated with the original KDA model and the solid ones with the 
modified model. 

As in Figure 3, the age markers are connected with dotted lines, and 
the member sources of each clan are indicated with the same symbols. Note 
that the ordinate of every source is its equipartition energy density derived 
from the observational data, while the tracks show the model cocoon energy 
density. This is an additional argument that the modified model better 
reproduce the observational data. In Figure 4, an uncertainty of the values 
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Figure 3: Evolutionary P—D tracks fitted for three clans of sources with 
evolving axial ratio AR (solid curves). The markers of the same predicted 
age on each curve are connected with dotted lines. The members of each 
clan are marked with different symbols as those in Figures 2a and 2b. Their 
actual age is indicated by a number behind the symbol. The dashed curves 
indicate relevant tracks but calculate with a constant AR, as in original KDA 
model. 

of both parameters is shown by the error bars. The error of u cq is given in 
Paper I (Table 2) and the error of E tot is estimated from the error of source 
volume (in Table 1 there). Again, the numbers mark the actual age of the 
corresponding source. It is worth emphasizing that the modified evolutionary 
u c -E tot tracks are steeper and curved in respect to those expected from the 
original model. This effect, related to the rate of adiabatic losses and inflation 
of the cocoon, is also discussed in Section 3.1. 

3 Discussion of the Results 

3.1 Dynamical Evolution of Individual Sources 

The data show that the source (cocoon) geometry, described by its axial 
ratio, depends on its estimated age. In Sections 2.3 and 2.4 we have pre- 
sented the evolutionary tracks of the 'fiducial' sources across the P-D and 
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Figure 4: Tracks of evolving energy fitted for the three clans. All the lines 
and symbols are the same as in Figure 3. 



■u e q-£tot planes, respectively, calculated with incorporation of the cocoon ax- 
ial ratio evolving with time into the analytical KDA model. A predicted 
evolution of each fiducial source, i.e. its tracks across the above planes, is 
verified by the age, luminosity, size, and the energy of a clan members (cf. 
Section 2.2). Therefore, we consider those sources as an individual source 
observed at different lifetime epochs, and the relevant fiducial source is its 
model representation. 

Figures 3 and 4 show that for each of the three clans our tracks are 
steeper than those resulting from the original KDA model, and much better 
corresponds to the observed parameters of the clan sources. Our P-D tracks 
strictly resemble those published by Blundell et al., as both arise from the 
dynamical models without assumption of the self-similar evolution of the 
source cocoon. Also the tracks of our clans across the u cq -E tot plane are 
steeper than those predicted by the original KDA model for sources with 
constant axial ratios. Moreover, the steepening increases throughout the 
source lifetime. This is provided with the time-dependent Phc ratio and 
non-constant inflation of the cocoon applied in our analysis. 

In the KDA model the synchrotron emitting particles flow from the jet 
head into the lobe (cocoon) whose pressures are in a constant ratio through- 
out the sources' lifetime. Thus the adiabatic losses in this model do not 
increase with time. In the BRW model the population of emitting particles 
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adiabatically leave the constant-pressure hotspot area into the cocoon whose 
pressure decreases throughout the sources' lifetime causing the losses due to 
its adiabatic expansion increase with time. It is interesting that we get sim- 
ilar tracks from a different approach, subsequently supporting predictions of 
their model. 

The observed distributions of the clan sources across the P D and u eq - 
E tot planes give further arguments after the necessity of accounting for the 
evolving axial ratio of the cocoon in studies of the dynamical evolution of 
individual radio sources, and allow to conclude that a departure from the self 
similar expansion in large and old sources is very probable. 

All the analyzed clans include giant sources, thus a dynamical evolution 
of radio sources from typical ('normal') size to 'giant' size is in accordance 
with predictions of the analytical KDA and BRW models. This supports our 
statistical findings that giants are old sources with high enough jet power 
evolved in relatively low-density environment. 



3.2 Deceleration/ Acceleration of the Cocoon Expan- 
sion Speed 

The expansion speed of the cocoon along the jet axis can be directly derived 
dividing length of the jet [equation (1) in Paper I] by time. Indicating g = 
Qo/po we have 

V h = L-Jt OC t(/*-2)/(5-/jy/(5-/J) 

In the above relation, taken directly from the KDA model, an acceleration 
of the cocoon expansion speed is admitted if /5 > 2. The value assumed 
in Paper I ,3=1.5 implies Vh oc t -a143 . However this relation is valid for 
the cocoon self-similar expansion only. As there is a doubt about validity 
of this assumption for the large and old radio sources, the expansion speed 
can rather be derived from the cocoon time-evolving volume of cylindrical 
geometry given by 

V c (t,g) = 2^[Lj(t,ff)] s oc ^7(5-/3)^3/(5-,) _ (1) 

(cf. equation (6) in Paper I). Inserting L } = D'/2 and i? T = AR/2 into Equa- 
tion (1) and taking 1/3 root of it, we have D(t,g) oc (ARf/H^ 1 ^-^ y/(s-£) . 
But the sample statistics indicate that AR also evolve with time and the 
ratio of Qo/po- Substituting 

AR(t : g)cxt - 16±om g°- 082±004 
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obtained from fitting a surface to the values of AR over the t—g plane and 
(5 = 1.5 into Equation (1), then dividing it by t, now we have 

V h (t, 9 ) OC t -0.04±0.04g0.34±0.03_ (2) 

The negative power exponent of t still implies the cocoon expansion-speed 
deceleration with time, though slower than that for the self-similarity ex- 
panding cocoon. 

Nevertheless, a slow acceleration effect is seen in the time evolution of 
the fiducial sources representing the clans of a few actual sample sources in 
which each source in the clan with almost identical values of Qo and po is 
observed at a different age (cf. Section 2.2 and Figure 2a). Below we analyse 
whether it can be real. 

Independent arguments after an accelerated expansion come from the 
well known hydrodynamical considerations. From Scheuer (1974), it has 
commonly been accepted that the ram pressure of the external gas behind 
the head of jet is 

p a V h 2 « Qo/(A h v iet ), (3) 

where p a is the external density, is the head working surface area, and 
the jet bulk velocity v- jCt is commonly assumed to be close to the light speed 
c. In this approach Vh is a function of the root square of the ratio Qo/p a and 
the reciprocal of the linear size (diameter) of the working area. Qo/Pa si m ~ 
ply represents an effectiveness of the jet propagation across the surrounding 
medium. We can consider this ratio at the core radius a and at the source 
leading head where it is Qo/po and Qo/Pend = Qo/po(-D/2ao) _/3 , respectively. 
Both ratios are derived from the model. Values of the first ratio for giant 
sources are not very different from those for other sample sources being lit- 
tle smaller than corresponding values for high-redshift sources, comparable 
with those for low-redshift sources, and evidently higher from those for low- 
luminosity ones. The second ratio, as a function of is the highest for 
giant sources. 

The area A^ at a given radius from the AGN centre can be determined 
from high-resolution observations of the hotspots. For example, VLBI obser- 
vations of selected 'compact symmetric objects' (CSO), possibly progenitors 
of classical double FRI and/or FRII sources give a linear size of the working 
area of about a few parsecs at a radius of about 50-hlOO pc (cf. Owsianik 
et al. 1998). Many VLA observations of hotspots in sources of size 10-i-20 
kpc show that their sizes are not larger than about 500 pc, while the hotspot 
sizes in the lobes of large sources do not exceed 15 kpc. Schoenmakers et al. 
(2000) in their study of giant radio sources have assumed a typical hotspot 
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size of 5 kpc at a radius of 0.5 ~ 1 Mpc. After the discovery of the third 
largest radio galaxy J1343+3758 (1343+379 in Tables 1 and 2), Machalski 
and Jamrozy (2000) determined a working surface area with a diameter of 
~13 kpc in one of its lobes. 

From the data in our sample sources we have (Qo/po) = 10 61 - 9±013 Wm 3 
kg -1 and (Qo/Pcnd) = 10 64 - 9±a19 Wm 3 kg -1 . Dividing the first ratio by 
Ah PS (2 ± 0.7) • 10 38 m 2 (the working area for a hotspot diameter of 0.5±0.1 
kpc) and the second ratio by Ah ~ (8 ± 6) ■ 10 40 m 2 (for a hotspot diameter 
of 10±4 kpc), we obtain Vh of 0.12 and 0.19c, respectively, which confirms a 
possible acceleration of the cocoon expansion speed. The above speeds, too 
high with respect to an upper limit of Vh (i.e. the quotient D/2t), can be 
caused by either an overestimation of the jet power or too low internal density 
of the core and of the external density at the end of jet. An adjustment of 
f3 will not work because some increase of j3 can compensate for too low 
Po, but will dramatically lower p en d. Since the values of A\ l adopted for 
this calculation may be overestimated rather than underestimated, only a 
decrease of Q can fit the observed expansion speeds. 

Indeed, the values of Qo derived from equation (3) base on the ram pres- 
sure considerations in the overpressured source model A of Scheuer (1974) 
and its further modifications (e.g. Begelman and Cioffi 1989; Loken et al. 
1992; Nath 1995). They all are self-similar models of the Carvaldo and O'Dea 
(2002) type I which describe the source dynamics only. If the source ener- 
getics and especially the energy losses are taken into account (the type III 
models; e.g. KDA, BRW), the significantly higher values of Q are implied. 
Therefore, the above calculation of Vh values must be based on Q Q values 
determined within the same dynamical model of type I. If this is the case, 
the controversy about too high expansion speeds will disappear. Using the 
Nath's (1995) model in which 100% of the twin-jets energy is converted into 
radiation (i.e. 2Q t = U cq ), we find that Qo values fitted with that model 
are about ten times lower than the KDA values. Therefore, the Vh values 
calculated as above but with the Q (NATH) values decrease to about 0.04 
and 0.06c, i.e. to the expansion velocities quite plausible. 

4 Conclusions 

The following conclusions can be drawn from the above analysis: 

(1) The sample sources, having similar model parameters Q and p but 
different age, size, axial ratio and luminosity (called a 'clan'), fit satisfactory 
the modified Kaiser et al.'s (1997) [KDA] model tracks on the P—D and u cq — 
E tot planes calculated for a 'fiducial' source being the model representation 
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of each clan. The derived tracks appear to be much steeper than those 
provided by the original KDA model; the best fit is achieved by allowing an 
evolution of the sources' cocoon axial ratio with time. The derived tracks 
are compatible with those expected from the more sophisticated Blundell et 
al.'s (1999) model. 

(2) The evolution of the energetics of sources on the u cq — E tot , predictable 
from the model, form another characteristic of their population which can 
be constrained by observational data. 

(3) In the case of clan sources, we found a slow acceleration of the average 
expansion speed of the cocoon along the jet axis. This effect is predicted with 
the model modified as above, and is caused by a systematic increase of the 
ratio between hotspot and cocoon pressures (Vh c ) with the age of source and 
decreasing density of the external environment according to the King's law 
model. The statistics suggest that this acceleration is also dependent on an 
efficiency of the jet propagation through the surrounding medium. 
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